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Abstract 

We point out that the anticommutation properties of the Dirac matrices can be derived with- 
out squaring the Dirac hamiltonian, that is, without any exphcit reference to the Klein-Gordon 
equation. We only require the Dirac equation to admit two linearly independent plane wave solu- 
tions with positive energy for all momenta. The necessity of negative energies as well as the trace 
and determinant properties of the Dirac matrices are also a direct consequence of this simple and 
minimal requirement. 
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I. INTRODUCTION 



Many textbooksi'^i'^i^'^i^ derive the Dirac equation for a free particle of mass m following 
the method used by Dirac himself in his 1928 paper-. This method involves two steps. 
First, one admits that the wave function should be a multi-component object, as in the 
non-relativistic theory of spin, and that its time evolution is ruled by a partial differential 
equation of first order in both the time and space derivatives. Working in a system of units 
where ^ = c = 1, we have 



where the Dirac hamiltonian Hd is defined by 



Hd = cx- (-iV) + j3m 

= E«^Ha^^) + /^-- (2) 
k=i ^-^ 



In this equation, the a^'s and /3 are constant hermitian matrices. In the second step, one 

'df 



'squares' Eq. ([T]) by acting on both sides of it with the operator ij-. This yields 



= Hl^- (3) 

Then one requires H^) to be identical to the operator — A + m^, thereby ensuring that each of 
the components of \l/ satisfies the Klein-Gordon equation. This implies the anticommutation 
relations 

ai(3 + [3ai = 0, 

f3' = 1. (4) 

Starting from these results, one usually proceeds by showing that such matrices indeed exist 
when \1/ is a four-component object and then one 'finds' that Eq. ([T]) admits both positive 
and negative energy plane wave solutions. So the Dirac equation does not solve the 'problem 
of negative energies' which appears when studying the Klein-Gordon equation. 

However, it is difficult to be immediately convinced that this fact is not a mere conse- 
quence of the requirement appearing in the second step of the above derivation. Thus, to 



discard any doubts about the necessity of negative energies, it would be more satisfactory 
to avoid the squaring of the hamiltonian Hd. In this paper, we show that this is indeed 
feasible. 

II. NECESSITY OF NEGATIVE ENERGIES 

In order to implement the program outlined at the end of the previous section, we require 
Eq. ([1]) to admit two linearly independent plane wave solutions of the form 

^(x,t) =M(p)e*(P-^-^p*\ (5) 

where E-p is the positive energy associated with a free particle of momentum p, that is, 

£;p = +(p^ + m2)V2. (6) 

Since our aim is to describe spin 1/2 particles such as electrons, this requirement is both 
natural and minimal. By inserting Eq. ([5]) into Eq. ([1]), we obtain 

Epn(p) = /iz)(p)m(p), (7) 

with 

^^(p) = • P + Z?"^- (8) 

Note that hnip) is a matrix of numbers whereas He, is a matrix of differential operators. 
Since we obviously discard the solution 'u(p) = 0, we see, from the above requirement, that 
Ep should be a double root of the eigenvalue equation pertaining to the matrix /iz)(p)- Thus, 
if we introduce the characteristic polynomial of hoiv) 

P^E) = dtm[E - hn{p)], (9) 

we should have 

Pn{Ep) = 0, (10) 

and 

^n(^p) = 0, (11) 

where is the derivative of P„ with respect to E. The index n in these equations stands 
for the degree of Pn{E) or, equivalently, for the number of components of the wave function 
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Let us now try to satisfy Eqs. (fTOj) and (fTT!) within a two-component theory (n = 2). We 
have 

P„(E) = P2(^) =E^ + ci(p)E + co(p), (12) 

where the coefficients Ci and Cq are polynomials homogeneous in m and the components pi, 
P2, Ps of the momentum p. Eq. (fTTj) yields 

2Ep + ci(p)=0. (13) 

It is not possible to satisfy this equation for all momenta since the square root Ep cannot 
be expressed as a polynomial. Thus, a two-component theory is immediately ruled out. So, 
let us try a Dirac equation with three components. Now, we have 

Pn{E) = P,{E) = E^ + C2{p)E^ + ci{p)E + co(p), (14) 

where the coefficients C2, Ci and Cq are again polynomials homogeneous in m and the com- 
ponents of the momentum p. Eqs. f llUj) and fllip yield 

El + C2ip)El + ci(p)Ep + co(p) = 0, (15) 

3El + 2c2ip)Ep + ciip) = 0. (16) 

Again using the fact that Ep cannot be expressed as a polynomial, we see that these equations 
imply 

El + c,{p) = 0, (17) 

C2(p)^J + Co(p)=0, (18) 

3El + c^{p)=0, (19) 
C2(p) = 0. (20) 

Eqs. (1171) and f|T9|) lead to E'p = for all momenta. This is not possible and, as a con- 
sequence, a three-component Dirac theory is also ruled out. Finally, let us turn to a four- 
component theory. Now, 

Pn{E) = P^{E) = E^ + cME^ + C2{p)E^ + ci{p)E + co(p), (21) 

where our notations are similar to those used above in the two- and three-component cases. 
Eqs. (dUD and ^ yield 

E^ + cMEl + C2{p)El + ci{p)Ep + co(p) = 0, (22) 
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4El + 3cs{p)El + 2c2(p)£;p + ci(p) = 0. (23) 

These equations imply 

E^ + C2{p)El + co{p)=0, (24) 



Csip)El + c^ip) = 0, (25) 
2Ej + C2(p)=0, (26) 



3cMEl + c^{p)=0. (27) 



From Eq. (!26l) . we obtain 



C2(p) = -2^;^. (28) 
Inserting this expression into Eq. ( l24l) yields 

co(p) = i^J. (29) 
Finally, comparing Eqs. (l25l) and (1271) leads to 

ci(p) = (30) 

and 

C3(p) = 0. (31) 

If we insert these results back into Eq. fl2T|) . we see that the eigenvalue equation for /i_d(p) 
reads 

{E ~ Ep)\E + E^y = 0. (32) 

This shows that the positive energy solutions to the Dirac equation will always be accompa- 
nied by solutions with negative energy. To prove that the approach adopted in this paper is 
self-contained, we still have to derive the anticommutation relations (jlj). This is performed 
in the next section. 



III. DERIVATION OF THE ANTICOMMUTATION RELATIONS 

We now show that Eqs. ([28]), ([29D, ([30]) and ([3l]) do indeed imply Eqs. ([!]). We remark 
that once we have replaced Ep by its expression ([6]), all of these equations require some 
polynomial homogeneous in m and the components of p to vanish identically, that is for 
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all momenta. This is possible only if all the polynomial coefficients are zero. We shall rely 
repeatedly on this remark in what follows. 
From Eqs. and (^li), we obtain 

C3(p) = -TrihDip)) 

3 

= Y.PkTr{ak)+mTr{P), (33) 

k=l 

where the symbol Tr denotes the trace. Thus, Eq. (l3T]) implies 

Tr(ai) = Tr(a2) = Tria^) = Tr{(3) = 0. (34) 
Eqs.© and (Ell) also yield 

Co(p) = dtm{hD{p))- (35) 

Inserting this expression into Eq. fl2^ and considering the terms in p\, pi, and leads 

to 

dtm{ai) = dtm{a2) = dtm{a^) = dtm{f3) = 1. (36) 

To make things simpler, it is convenient to work in a representation where the matrix 
P is diagonal. Note that Eqs.llMl) and ( l36i) are representation independent. Consider the 
terms in in Eq. fl30|) and in in Eq. fl28l) . They yield 

Al/322/333 + Pll/322/3u + /3ll/333/?44 + (322(333/3u = (37) 

and 

+ Al/533 + /3n/3u + /322(333 + (322pU + /?33/?44 = "2, (38) 

respectively. Combining these equations with 

dtm{(3) = /3n/922/333/944 = 1, (39) 

(see Eq.(l36!)). we obtain 

(1 + Ai)(l + /522)(1 + /333)(1 + M = (40) 

and 

(1 - /5n)(l - /322)(1 - PssKl - I3u) = 0. (41) 

These equations show that one of the eigenvalues of (3 is equal to +1 and another to — 1. 
Let us assume that Pu = +1 and P33 = —1. Taking Eqs. flMl) and (|36|1 into account, this 
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implies = —P22 and P22 = ±1. We shall assume that P22 = +1 and thus ^4^4 = — 1. We 
do not have to consider other choices for the diagonal elements to be put equal to +1 or —1 
since this would correspond to a mere rearrangement of the lines and columns of /3. Thus, 
we have 

^10 
10 
0-10 
-1 



P 



\ 



\ 



(42) 



Obviously, 



and it is easy to show that this equation implies 



(43) 



(44) 



1 = 1,2,3). 

Indeed, let us just imagine that we perform, on all the Dirac matrices, a unitary trans- 
formation which brings ai, say, into diagonal form. We expect that the matrix /5 will no 
longer be diagonal but Eq. (H^ will remain true because it is representation independent. 
We now proceed for ai as we did above for /?, that is, we concentate on the terms in 
in Eq. fl30|) and in pf in Eq. fl28l) . This will lead us to = 1 which is also representation 
independent. Proceeding in this way for 0:2 and as, we prove the other identities in Eq. (l44l) . 
This trick can be used each time we establish a representation independent identity even if 
we arrived at that identity within a particular representation. In what follows, we go back 
to the representation in which /5 is given by Eq. fl42l) and we derive the structure of ai in 
that representation. For the moment, we drop the index 1 to simplify our notations. Thus, 
a stands for ai. Consider Eq. (1281) . The terms in mpi and in pf give 

— an — a22 + 033 + "44 = (45) 

and 



«11«22 + «11«33 + ailCXu + «22«33 + a22«44 + «33«44 
I |2 I |2 I |2 I |2 I |2 I |2 c 

— |ai2| — laisl — |ai4| — |a23| — |«24| — |«34| = —2 
respectively. Consider now the terms in m'^pl in Eq. (l29l) . they give 

«lia22 — «lia33 — aiia44 — a22a33 — a22«44 + «33«44 



(46) 



I |2 I |2 I |2 I |2 I |2 I |2 

— |ai2| + [aisl + |ai4| + |a23| + |tt24| — F34I = 2. 



(47) 



Adding Eqs.(ll6]) and dH]) yields 



aiia22 + a33«44 — |ai2|^ — |a34p = 0. 



(48) 



On the other hand, comparing Eg. (1451) with Eq. (lM|) yields 

^22 = -an 

and 

0:44 = -033. 

If we insert these results back into Eq. fHHl) . we obtain 



«ii + (^33 + I"i2p + |a34p = 0. 



Thus, we have 



an = "22 = "33 = "44 = "12 = "34 = 0, 

and, restoring the index 1, we see that the matrix ai has the following structure: 



(49) 



(50) 



(51) 



(52) 



\ 



(ai)i3 (ai)i4 
(ai)23 (ai)24 

(«l)23 

\(ai)t4 ("1)24 y 
where the non-vanishing elements are restricted by the condition 



(53) 



|(ai)i3|' + |(ai)i4|' + |(tti)23|' + |(«i) 
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2. 



(54) 



Actually, this equation tells us nothing new since it can be derived from Eq. flHl) . An 
analogous proof shows that the matrices 02 and have also this structure. We note that 
the structure of the a^'s and of P (see Eq. fH2|) ) imply 



aip + pai = (i = 1,2,3), 



(55) 



as can be checked simply by performing matrix multiplications. Since these equations are 
representation independent, we conclude, using the trick described after Eq. fH^ . that we 
should also require 

aiaj + ajai = (i 7^ j = 1, 2, 3). (56) 
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Eqs.( l43p . (jBD, (155|) and (156|1 are the anticommutation relations we were looking for. It is 
easy to check that Eqs. fl28|) . fl29|) and fl30l) do not give rise to additional restrictions on the 
Dirac matrices. As an example, consider the terms in piP2 in Eg. (1281) . They impose 



and, indeed, requires nothing new. 

IV. SUMMARY AND COMMENTS 

In this paper, we have provided a method to derive the anticommutation properties of 
the Dirac matrices without relying on the squaring of the Dirac hamiltonian. We have 
only required Eq.([T]) to admit two linearly independent plane wave solutions with positive 
energy for all momenta. At an early stage in the derivation, we have seen that, despite 
this conservative requirement, it was not possible to rule out negative energy solutions, 
thereby establishing that these are not an artefact of the standard derivation. It might 
also be interesting to note that, within the method described in this paper, the trace and 
determinant properties of the Dirac matrices appear in the course of the derivation and not 
as by-products of the anticommutation relations. Finally, a few comments are appropriate 
concerning our proof of the impossibility of a two-component Dirac equation. As is well 
known, such an equation appears in a space-time with less than three space dimensions^ or 
in the study of massless fermions, where it is known as the Weyl equation^i^. This by no 
means contradicts our assertions. Indeed, no impossibility arises in a two-component theory 
if one only requires the Dirac equation to admit a single plane wave solution with positive 
energy for all momenta. In that case, only Eq. (fTOl) with n = 2 has to be imposed and this 
yields 



(ai)i3("2)i3 + (ai)i3(a2)l3 + (ai)i4("2)i4 + (ai)i4(a2)M 

+ («l)23("2)23 + (ttl)23(«2)23 + (ai)24(«2)24 + (ai)24(«2)24 = 0- 



(57) 



This equation can be written as 




(58) 




(59) 



This equation implies 




(60) 
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and 



ci(p) = 0. 



(61) 



From these equations, we see that the eigenvalue equation for /id(p) reads 



{E-Ep){E + Ep) = 0. 



(62) 



Thus, we have a plane wave solution with positive energy and another with negative energy. 
As a consequence, in a two- component theory, the 'twofold degeneracy' only corresponds to 
the existence of antiparticles. The derivation of the properties of the Dirac matrices (actually, 
of the Pauli matrices, since we are now in a two-component theory) can be performed as in 
the previous section and will not be repeated here. 
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